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Abstract. We briefly review the recent results obtained by Fermilab experiment E791 on the Dalitz 
plot analysis of the decay D+ K^k^k^, where indication for a light Kn scalar resonance, the fc, 
was found. We also present preliminary studies providing further information on the phase behavior 
of the scalar components at low mass, supporting the previous indication for the K. 

INTRODUCTION 

Nowadays, with the rising statistics of chann samples, the decays of charm mesons can 
be seen as a new source for the study of light meson spectroscopy, complementary to 
that from scattering experiments, and can be particularly relevant to the understanding 
of the scalar sector. 

Recently, Fermilab experiment E791 reported on the indication of a light scalar Kn 
resonance ||l|], k, based on a full Dalitz plot analysis of the decay K^n^n^. The 

measured Breit-Wigner mass and width of this state were found to be 797 ± 19 ±43 
MeV/c^ and 410 ±43 ±87 MeV/c^, respectively. Fermilab E791 has also observed 
evidence for the scalar o meson in D+ n^n^K^ decays L2J and measured /o masses 
and widths in n^n^n^ decays 

Here, we briefly review the K^Tt^Tf^ methodology and the results obtained. 

Moreover, we present new, preliminary studies providing further information on the 
phase variation of the scalar sector at low mass. We first present a study of the interfer- 
ence of ^*(890) and underlying amplitudes as an attempt to measure the phase motion 
of the scalar low-mass Kit amplitude. A more extensive study comes afterwards, where 
we try a set of fit models inspired by the LASS parametrization for the S-wave Kit am- 
plitude [4J. The main idea is to check whether a unitary, isolated, S-wave Kit amplitude 
can fit the Z)+ K^it^it^ Dalitz plot, together with the other (higher-spin) intermedi- 
ate states. We show that this approach is not adequate for the data, and present similar 
parametrizations where final- state interactions are included for the S-wave components, 
as is usually done in Dalitz plot analyses. In particular, we present a result offering fur- 
ther support for a rapidly- varying phase at low Kk mass, consistent with the previous 
indication of a light, broad K in D+ K^it^it^ . 



REVIEW OF THE D+ K' 71+71+ RESULTS 



From the original 2 x 10^*^ events collected by Fermilab E791 jsl], and after reconstruc- 
tion and selection criteria, we obtained the K^Tt^ 7t^ mass spectrum shown in Figda). 
The sample of 15,090 events shown in the crosshatched area is used for the Dalitz plot 
analysis of K^n^n^, from which about 6% is due to background. The Dalitz 

plot of these events are shown in Fig[Ttb), where the axes represent the two invariant 
mass squared combinations Si2 and 513 (the kaon is labeled particle 1, the two pions are 
particles 2 and 3). 

For the Dalitz plot analysis, an unbinned maximum-likelihood fit was performed with 
probability distribution funtions (PDF's) for both signal and background. The signal 
PDF was written as the square of the total physical amplitude £/ and it is weighted by 
the acceptance across the Dalitz plot and by the level of signal to background for each 
event. 

Many K* resonances could contribute to the final state by forming quasi 2-body 
intermediate states K*7i:^, as well as a possible non-resonant (NR) 3-body decay. Due 
to the possibility of rescattering of these intermediate states from final state interactions 
(FSI), a general ansatz, widely use in Dalitz analyses of D decays to describe the total 
decay amplitude, is given by 

^ = aNR^'^'^M^R + £ an^'^"< ■ (1) 

M=l 

This equation uses Lorentz invariant amplitudes £/„ to describe the individual resonant 
and NR processes, with relative strenghts given by the magnitudes a„ and strong phases 
dn to account for FSI. The NR amplitude ^^^nr is usually assumed constant and set to 

1. The resonant amplitudes are written as = BWn F^^ Fn^ where 5W„ is 

the relativistic Breit-Wigner propagator, 5W„ = \n^^ — irp- — imQY{m)\ \ the quantities 
Fd and Fr are the Blatt-Weisskopf damping factors respectively for the D and the Kn 

resonances and ^n'^ describes the angular function according to the spin J of the 
resonance. See JII]. The amplitudes are generally functions of both Dalitz variables 
512 and s\2,. Here, with two identical pions in the final state, each amplitude is Rose 
symmetrized £^ = i2^j[(12)3] +^,[(13)2]. 

The first approach to fit the data includes the known K* resonances plus a con- 
stant NR term. This is called Model A. It finds contributions from the channels: NR, 
with a decay fraction over 90%, followed by ^*(1430);r+, ^*(892);r+, ^*(1680);r+ 
and ^|(1430);r+; these values are in accordance with previous results from Fermilab 
E691[6] and Fermilab E687[7]. See fit fractions and relative phases in Table [T] There 
is an important destructive interference pattern, since all fractions add up to 140%. Im- 
portant disagreements appear between the fit to Model A and data, with = 2.7 (v 
being the number of degrees of freedom) in a binned version of the Dalitz plot. The 
discrepancies are found mainly at low Kn mass squared (below 0.8 GeV^/c^) and near 
2.5 GeVVc"^ lH. Thus, a model with the known Kn resonances, plus a constant NR 
amplitude, is not able to describe the D+ — > K^K^it^ Dalitz plot satisfactorily. 



(a) (b) 

FIGURE 1. (a) The Knn invariant mass spectrum. The filled area is background; (b) Dalitz plot 
corresponding to the events in the crosshatched area of (a). 



TABLE 1. Results without K (Model A) and with K (Model C). 



Decay 
Mode 


Model A: No K 
Fraction (%) Phase 


Model C: With K 
Fraction (%) Phase 


NR 

K*{%92)n+ 
^o*(1430)7r+ 
^:|(1430)7r+ 
^:*(1680)7r+ 


90.9 ±2.6 0° (fixed) 

13.8±0.5 (54±2)° 
30.6 ±1.6 (54 ±2)° 
0.4±0.1 (33±8)° 
3.2±0.3 (66±3)° 


13.0±5.8±4.4 (-11±14±8)° 
47.8±12.1±5.3 (187±8±18)° 
12.3 ±1.0 ±0.9 0° (fixed) 
12.5 ±1.4 ±0.5 (48 ±7 ±10)° 
0.5±0.1±0.2 (-54±8±7)° 
2.5 ±0.7 ±0.3 (28 ±13 ±15)° 



A second model (Model B) allows the mass and width of the scalar Kq (1430) to float, 
and includes Gaussian form-factors [1]. The resulting fractions and phases are similar to 
those of Model A (within errors) and the mass and width of (1430) are found to be 
consistent with PDG values. The fit improves but it is still unsatisfactory. 

A third fit model, Model C, is constructed by the inclusion of an extra scalar state, 
with unconstrained mass and width. The mass and width of the (1430) are also free 
parameters, and Gaussian form-factors are used as in Model B. Using this model, the 
values obtained are 797 ± 19 ± 43 MeV/c^ for the mass and 410 ± 43 ± 87 MeV/c^ for 
the width of the new scalar state, here referred to as the K. The values of mass and width 
obtained for the K^{U3Q) are respectively 1459 ±7 ± 5 MeV/c^ and 175 ± 12 ± 12 
MeV/c^, appearing heavier and narrower than presented by the PDG. The fit fractions 
and phases from Models C are shown in Table [T] Compared to the results of Model A 
(without k), the NR mode drops from over 90% to 13%. The kk^ state is now dominant, 
with a decay fraction about 50%. Moreover, the fit quality of Model C is substantially 
superior to that of Model A; the / v is now 0.73. 



From the results of the Dalitz plot analysis of — > K^n^n^ , one concludes that a 
conventional approach, including known Kk resonances and a constant NR term, cannot 
describe the data. The presence of a broad, light scalar Kn state, k, provides a very good 
description of the data and becomes the major contribution to the decay. Many other 
possibilities were tried with no sucess, for instance a toy model for which the extra state 
is represented by a Breit-Wigner without phase variation, vector and tensor hypotheses, 
different parametrizations for the NR channel, etc. 



FURTHER STUDIES 

As described above, we obtain a very good description of the D+ K^K^it^ data 
by the inclusion of the K in the decay amplitude. A scalar resonant phase behavior 
seems to be required at low mass. Nevertheless, it is desirable to investigate further, 
and independently of the Breit-Wigner hypothesis. 

Here, we present new studies in order to provide further information on the phase 
behavior of the scalar amplitude at low mass. A first study refers to the asymmetry of 
the ^*(890) and we show that there is no sensivity to the K phase motion at the ^*(890) 
mass in the crossed channel, although an indication of it can be seen at lower and higher 
mass. A second, more important study, refers to a comparison with the results of the S- 
wave Kk amplitude from LASS, composed of a non-resonant background (with effective 
range parametrization) and the ^q(1430). The main results we obtain are the inability 
to describe the data by imposing a Kit elastic unitarity constraint but, on the other hand, 
the features of the phase motion obtained after this imposition is released. 



/<:*(890) Asymmetry Study 

A suggestion has been made to us |8] to see whether the phase motion of the K 
could be inferred through the interference of the ^*(890) with the broad, underlying 
amplitude. From FigUJa) a strong asymmetry in the interference of the ^*(890) is 
evident in the Dalitz plot. This asymmetry is due to the interference with the dominant 
scalar component at low mass, i.e., the NR component in Model A (without k) or the 
broad K in Model C. In the following, we will explore this possibility. 

If one plots the phase of the k Breit-Wigner {5bw) with central mass and width 
obtained by E791 (797 MeV/c^ and 410 MeV/c^, respectively) one sees that there is 
a sharp increasing of the phase from threshold up to about 1 GeV/c^ (see Figl2|l. Above 
1 GeV/c^, the increase is slow. Now, we want to study the K phase motion in mi3 by 
inspection of the cos 9y distribution of ^*(890) {9y being the decay angle of the K in 
the center-of-mass frame of ^*(890) with respect to the D direction) within 0.8 and 
1.0 GeV/c^ in the crossed channel, m\2. By observing the Dalitz Plot, we see that for 
0.8 < mi2 < 1 .0 GeV/c^, the allowed kinematical region of the crossed channel is mi3 
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FIGURE 2. The Breit-Wigner phase 5bw for a scalar Kn resonance with mass 797 MeV/c^ and width 
410MeV/c2 

1 GeV/c^ Thus, we can anticipate that the study of the ^*(890) asymmetry will miss 
the main K phase variation in the crossed channel. 

From the data, this fact becomes evident. In Fig|3l we plot the cos By distributions 
(a) within the ^*(890) mass region and (b) for low and high Kk masses (outside the 
K*{%90) region). For all plots, we compare the data (points with error bars) to the result 
of our fit models with and without K. From FigOta) it is clear that it is not possible to 
distinguish between the models. Nevertheless, for low and high Kn masses shown in 
FigEtb) we see that the model without K gives a poor description of the distributions, 
while the model with k describes them adequately. 

As a complementary view, the asymmetry distribution is obtained by plotting cos Oy 
weighted by the number of events in each Kn mass bin between 0.8 and 1.0 GeV/c^. 
For E791 data, the distribution is shown in Fig|4|with error bars. In the same figure, both 
models with (Model A) and without K (Model C) are compared to the data distribution. 
Both reproduce very well the asymmetry observed for the ^* (890) . Thus, using only this 
information one would not claim the necessity of the K in the K^n^n^ decay. 

On the other hand, the effects of the K can be seen below the ^*(890) mass or further 
above it, as shown before. 

It is worthwhile to comment on the asymmetry reported by Fermilab FOCUS [9] in 
their analysis of the decay ^ K^n^jl^v^, where the ^*(890) was found to interfere 
with a small scalar component with a phase of 45°. There, no FSI phases are expected, 
thus a comparison with the phases of D+ K^n^n^, as suggested by Jl], could be far 
from being direct (see more discussion in the next section). 



It is very important to call attention here to one feature of the D+ K n'^n'^ decay. When looking at 
a low-mass region in one channel, let's say sn < 0.7 GeV/c^^ (just below K*{^90) mass) we are at the 
same time observing a higher region of sn, above 1.3 GeV^/c^. This is a result of having two identical 
pions in the decay. Thus, the presence of the K in this decay, especially its phase variation, would reflect 
itself at low and high Kn mass (respectively below 0.8 GeV/c^ and above 1.5 GeV/c^). 



With K M(K7t)<0.8 Without k 
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FIGURE 3. The cos By distribution: (a) within the K* (890) region (0.8-0.9 GeV/c^); (b) for low (< 0.8 
GeV/c^) and high (> 1.5 GeV/c^) Kn mass 
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FIGURE 4. The asymmetry distribution in cos 9v for Kn mass between 0.8 and 1 .0 GeV/c^. Error bars 
represent the data, Model A is shown as the dotted histogram, and Model C as the solid-line histogram. 

Comparisons with the LASS Kti Amplitude 

A set of studies was performed to compare the features of the S-wave Kn components 
in the — > K^n^n^ decays to the S-wave measured by LASS experiment [4] on 
elastic Kn scattering. It is interesting to discuss here the differences between these two 
environments. We will be comparing production to a scattering process. For D decays, 




FIGURE 5. The S-wave phase motion obtained by the LASS experiment (solid Une). The behavior of 
5nr is shown through the dashed Une. 



it is known that FSI can play an important role; explicitly there is the possibility of 
rescattering between quasi two-body states, RiTt^ R2^^ (where R12 are any of the 
Kn resonances). Thus the presence of the bachelor pion can have an influence on the 
process, in this case not being isolated from Kn. One consequence of this fact would be 
the loss of unitarity in the Kn system. This is what we want to check. 

The idea here is to try the same parametrization used by LASS to see whether the 
S-wave phase motion observed by them could represent the K^K^it^ data. Their 

S-wave (up to the inelastic threshold) had two terms: a non-resonant term, described 
through an effective-range form, and a relativistic Breit-Wigner for ^q(1430). No low- 
mass scalar resonance was observed. We will try this same kind of approach in two 
ways. First, we impose the exact form they used, which implies the imposition of Kn 
elastic unitarity. Then, we release this restriction. 

The amplitude of the S-wave used by LASS, translated to a D decay (instead of 
scattering) needs a phase space correction given by m/p* (mass of the resonance and 
decay momentum in the resonance frame, respectively). In the first fit. Fit 1, the relative 
strength and phase for NR and Kq are fixed by imposing the unitarity constraint on the 
Kk system, so the total S-wave amplitude is written as: 

= ^ sin 5s e''^ + e'^'^ , (fPn^ ^^^^^^ 
Pn m5 - mf2 - imor(mi2) 

with r(m) = (mo/m) (p* / p*'^)rQ and the NR term has the effective range form: 

cot dB = -^ + lrbp* (3) 
ap* 2 



The values obtained by LASS in their fit were 



a= 1.95±0.09c/GeV , = 1.76 ±0.36 c/GeV 
MK*{im = 1435 ± 5 MeV/c2 , r^*(i43o) = 279 ± 6MeV/c2 (4) 

From these values, we plot the S-wave phase motion they found from their data in Fig|5l 
The total amplitude describing the D+ — > K^k^k^ decay is written as a coherent 
sum of the above S-wave term (EqO and the other resonant states with spin > 1 . The 
relative strenghts and phases {a' s and 5' s below) are let free in the fit: 

j>\ 

^TOT = ase'^'s^s + £ a^e'^''^ (5) 

n 

When we fit the data with this model, we obtain a very bad representation. The /v is 
10.6. This is evidence for the fact that the imposition of the Ktz unitarity constraint for 
the S-wave does not work here: the S-wave Kk system is not isolated. 

Let us now drop the constraint of elastic unitarity in the Kit system, and use a 
model closer to what is the usual Dalitz plot approach, where the presence of FSI is 
included through extra strong phases for all contributing channels. As before, the scalar 
components included are the NR (parametrized as effective range) and the Kq{\A30) 
but each of them now has an independent magnitude and phase. Thus, the total decay 
amplitude for D+ K^k^tz^ is being written as: 

s^TOT = ^NRe'^^ ^ sin dB e'^^ + a^e^^S bWk* {mn) + (2^3) 
Pn 

j>i 

+ Y,ane'^"£^n. (6) 

n 

Using this approach, we try two fits: in Fit 2, we fix the parameters of the scalar 
components to the values obtained by LASS, shown in Eq|4l For Fit 3, these parameters 
(a, b, M^*(i43o), r^*(i430)) are allowed to float. 

Our results for Fit 2 are much better than those for Fit 1. Nevertheless they are very 
much like the results for Model A - the fit with a constant NR term, and without K. The 
quality of Fit 2 is still not good (X^/ v = 2.3) and the NR term appears with an enormous 
fraction of about 120%. 

Our results for Fit 3, on the other hand, are very interesting. We use the amplitude 
as in EqI3 but releasing the parameters a and b, and the mass and width of ^q(1430). 
In this case, we obtain a good description of the data. The values obtained for these 
parameters are: 

a = 4.58±0.32c/GeV , = -2.94 ±0.43 c/GeV 
^/fo*(i430) = 1459±6MeV/c2 , r^.(i43o) = 156± llMeV/c^ (7) 

Comparing to the best model, where K is included (Model C), we observe: the fit 
qualities are comparable (the model with K is slightly better); the mass and the width of 
the (1430) agree; the fraction of the non-resonant in Fit 3 (66 ± 7%) is similar to the 
sum of the NR and K in Model C; results for the other states agree. 
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FIGURE 6. NR phase in Fit 3 compared to Model C and to LASS 



A very interesting feature of Fit 3 is the negative value for b. In this case, the 
"effective-range" would be indicating a resonant behavior. To better show this effect, 
we plot in Figl^lthe 5b phase motion of the "NR" decay in Fit 3 for the fitted values of a 
and b. Note that the phase increases very rapidly from threshold. This effect is similar to 
the effect we found for K". Here in Fit 3, the "effective-range" phase would be replacing 
that for a constant NR and the K in Model C (also shown in the figure); the two curves 
are very similar, even more if one realizes that one is mainly the result of a Breit-Wigner 
shape for the K and the other comes from what would be an effective-range form (but 
turned out to have an impressive phase variation). For comparison, the curve of ds with 
the parameters obtained by LASS are also shown: an almost constant phase for the NR 
was observed. 

We see that the scalar sector at low Kn mass in K^n^n^ needs a rapidly- 

varying phase, consistent with the previous indication of the K state from UJ ■ 

CONCLUSIONS 

Using the Fermilab E791 data sample of Z)+ K^tz^tz^ , we present new, preliminary 
studies offering further information concerning the phase motion behavior of the scalar 
amplitude at low Kk mass. In particular, we showed that the S-wave phase motion 
as measured by LASS in a scattering process is not able to describe the E791 decay 
data. This result shows the significance of final-state interactions in D decays preventing 
a direct comparison of the phases in scattering and decay processes. Furthermore, an 
independent parametrization (without using a Breit-Wigner) provided a phase behavior 
consistent to what we found previously for the scalar low-mass K% amplitude as a K" 
meson. 



REFERENCES 



1. E. M. Aitala et al. (E791 collaboration), Phys. Rev. Lett. 89, 121801 (2002). 

2. E.M. Aitala et al. (E791 collaboration), Phys. Rev. Lett. 86, 770 (2001). 

3. E.M. Aitala et al. (E791 collaboration), Phys. Rev. Lett. 86, 765 (2001). 

4. D. Aston et al. (LASS Collaboration), Nucl. Phys. B296, 493 (1988). 

5. S. Amato et al., Nucl. Instrum. Meth. A324, 535 (1993); S. Bracker and S. Hansen, |h^-ex/02 1 0034{ 
S. Hansen et al., IEEE Trans. Nucl. Sci. 34, 1003 (1987) 

6. J.C. Anjos et al. (E691 Collaboration) , Phys. Rev. D 48, 56 (1993). 

7. PL. Frabetti et al. (E687 Collaboration), Phys. Lett. B 331, 217 (1994). 

8. W. Ochs, private communication. 

9. J. M. Link et al. (FOCUS Collaboration), Phys. Lett. B 535 43, (2002). 

10. W. Dunwoodie for the LASS Collaboration, private communication. 



